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Abstract: With spintronics applications in mind, we use irreversible thermodynamics to derive the rates of entropy 
production and heating near an interface when heat current, electric current, and spin current cross it. Associated with 
these currents are apparent discontinuities in temperature (AT), electrochemical potential (A/i), and spin-dependent 
"magnetoelectrochemical potential" (Afif^i). This work applies to magnetic semiconductors and insulators as well as 
metals, due to the inclusion of the chemical potential p, which usually is neglected in works on interfacial thermodynamic 
transport. We also discuss the (non-obvious) distinction between entropy production and heat production. Heat current 
and electric current are conserved, but spin current is not, so it necessitates a somewhat different treatment. At low 
temperatures or for large differences in material properties, the surface heating rate dominates the bulk heating rate near 
the surface. We also consider the case, noted by Rashba, where bulk spin currents occur in equilibrium. Although a 
surface spin current (in A/m^) should yield about the same rate of heating as an equal surface electric current, production 
of such a spin current requires a relatively large "magnetization potential" difference across the interface. 

PACS Nos.: 05.70.Ln,05.70.Np,67.40.Pm,73.40.Cg 

Resume: Avec applications dans 1' esprit de spintronics, nous employons la thermodynamique irreversible a obtenir 
les taux de production d'entropie et de chauffage a proximite d'une interface lorsque la chaleur actuelle, le courant 
electrique, et courant de spin la traverser. Associes a ces courants sont discontinuites apparentes de la temperature (AT), 
potentiel electrochimique (A/i), et dependant du spin potentiel "magnetoelectrochimique" {Ap.-f,i). Ce travail s'applique 
a semi-conducteurs magnetiques et isolants ainsi que des metaux, due a I'inclusion de la potentiel chimique /i, ce qui 
est generalement negligee dans les travaux sur les transports thermodynamique interfaciale. Nous discutons egalement 
de la distinction (non evidente) entre la production d'entropie et la production de chaleur. Chaleur actuelle et le courant 
electrique sont conserves, mais n'est pas courant de spin, il necessite un traitement quelque peu different. A basse 
temperature, ou pour de grandes differences dans les proprietes du materiau, la vitesse de chauffage de surface domine la 
Vitesse de chauffage en vrac pres de la surface. Nous considerons egalement le cas, a note par Rashba, ou les courants de 
spin en vrac se produire a I'equilibre. Meme si un courant de spin de surface (en A/m^) devrait donner environ le meme 
taux de chauffage d'une surface egale de courant electrique, la production d'un tel courant de spin necessite un potentiel 
relativement important "aimantation difference" entre 1' interface. 



1. Introduction 

It is well-known that apparent voltage and temperature dis- 
continuities, determined by extrapolation from the bulk, ap- 
pear at interfaces in the presence of heat or electric current. 
For small currents, these discontinuities are proportional to 
the heat or electric current. For heat current, the coefficient of 
proportionality is known as the thermal boundary resistance, 
and was first studied at low temperatures by Kapitza for the 
solid-liquid ^He interface [1^]. For electric current, the co- 
efficient of proportionality is known as the surface resistance, 
or specific resistance [5]. In principle, there can also be off- 
diagonal terms, corresponding to a discontinuity in the tem- 
perature causing an electric current [6]. There also are more 
recently, spin-dependent conduction effects across surfaces, as 
studied, for example, in [6-8]. 

Johnson and Silsbee [6] studied the surface and bulk trans- 
port coefficients for these phenomena, but without considering 
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details of the non-conservation of the spin current due to spin- 
flip processes, and did not study the rate of heating near the sur- 
face. The present work considers these non-conservation phe- 
nomena, which require more refined considerations than when 
they are not present. (The terms "magnetoelectrochemical po- 
tential" and "magnetization potential" that appear in the theory 
were first employed in [6]; they are made more precise below.) 
The present work also includes the chemical potential /i of 
the charge-carriers, neglected by [6], which considered metals. 
The present results are more general; /i is negligible (compared 
to electrical potential energy) for metals but not necessarily for 
semiconductors or insulators, where small changes in carrier 
density can have a large effect on fi. 

The present work also considers the conditions under which 
bulk heating dominates surface heating, and vice-versa. It also 
explicitly considers the near-surface region. Note that heat- 
ing implies entropy production, but not the converse; Sect. IIB 
presents some considerations on this matter. 

For specificity, note that if there is an extrapolated tempera- 
ture discontinuity AT across a solid-solid interface, the energy 
flux (and heat flux) are given by 

\f\^hK\AT\, (1) 
where Hk > is the Kapitza, or thermal boundary, conduc- 
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tance and Rk = is the Kapitza, or thermal boundary, re- 
sistance. Normally this expression is simply written down, [6] 
being an exception (with a very different notation). (As usual, 
heat flows from hot to cold, and we assume no voltage drop 
across the interface.) Entropy is produced but no heating oc- 
curs, as thermal equilibration occurs by heat flow. 

In addition, when there is an extrapolated voltage discon- 
tinuity AF across a solid-solid interface, the electric current 
flux is given by 



\3\=9\^V\, 



(2) 



where g is the surface electrical conductance. Likewise, this 
equation normally is simply written down, although we note 
[6-8]. To our knowledge, in none of these cases is the relation- 
ship to heating made. (As usual, current flows from high to low 
voltage, and we assume no temperature drop across the inter- 
face.) In this case there is both heating and entropy production. 

Although both heat current and electric current are associ- 
ated with a rate of entropy production Ss, only the latter is as- 
sociated with the rate of heat production. Specifically, the rate 
of heating (per unit area) associated with an interface is given 
only by 



Hi = |j(AF)|, 



(3) 



a simple generalization of the bulk rate of heating (per unit 
volume) of j • E. Furthermore, when there also is a spin current 
J" (in units of lis), there is an additional rate of heating (per 
unit area) associated with an interface, given below, as the third 
termin eq.(50). 

The fundamental principle that ensures (3) is the non-decreasinj 
nature of the entropy of the overall system, just as that same 
principle ensures the corresponding relation (2) between elec- 
tric current and voltage gradient; both yield current flow from 
high to low voltage. 

General references on irreversible thermodynamics are [9, 
10], and appropriate sections in the thermodynamics texts by 
Callen and by Morse [II, 12]. A number of more recent ap- 
proaches and applications, closer to the spirit of the present 
work, are also available [13-20]. The present work does not 
consider length scales so small that ordinary heat conduction 
(e.g., the Fourier law) is not expected to hold, due either to 
classical or quantum size effects [21]. Although very general, 
irreversible thermodynamics only applies to systems near equi- 
librium; it does not treat situations Uke the free expansion of a 
gas. 

The present work does not consider systems for which the 
interface is a distinct thermodynamic system. For a system 
that is out of thermodynamic equilibrium a surface temperature 
may not be a well-defined quantity. (See [22, 23] for molecular 
dynamics simulations of heat flow across an interface, which 
show a sharp temperature jump at the atomic level.) Moreover, 
thermometers that measure different properties, but are cali- 
brated in the bulk, need not read equivalent temperatures near 
the surface. This is because near surfaces the thermal distribu- 
tion function is not defined solely in terms of thermodynamic 
properties, but also in terms of surface solutions of the trans- 
port equation [24, 25]. 

Section 2 considers heating in a single charge-carrier sys- 
tem (or even a system with no charge-carriers) and obtains, in 



addition to the well-known bulk heating rate, the surface heat- 
ing rate. Section 3 considers a two charge-carrier system (with 
spins in mind), obtaining the bulk heating rate (including equi- 
librium spin currents, which, if non-uniform, must have zero 
divergence) and the surface heating rate. Section 4 provides a 
brief summary. 

A number of additional considerations are given in the Ap- 
pendices. The approach taken in the present work employs in- 
tensive thermodynamic variables associated with the energy- 
maximum principle (e.g., T, P, ji), and considers the case 
where certain of the extensive thermodynamic variables are not 
conserved. On the other hand, the approach taken in [11, 12] 
employs variables associated with the entropy-minimum prin- 
ciple (e.g., I/T, —P/T, — /i/T), but considers only the case 
where all of the extensive thermodynamic variables are con- 
served. Appendix A shows how to use the entropy-minimum 
variables to obtain the irreversible thermodynamics when the 
extensive thermodynamic variables are either conserved or not 
conserved. Appendix B shows that, in the long- wavelength Umit, 
the entropy of a sound wave is zero. Appendix C compares 
the notation of the present work with that of [6]. Note that 
Ref. [6] considers metals, whereas the present work contains 
results that are also applicable to magnetic semiconductors. 

The same approach can be applied to interfaces involving 
electrons and holes, rather than up-spin and down-spin elec- 
trons. Just as the spin current is not conserved, because up- 
spins and down-spins can flip, so too the difference between 
the electric currents due to electrons and holes is not conserved, 
due to electron-hole recombination. 



2. Single Carrier Systems 

The spirit of irreversible thermodynamics is to write down 
the differential of the energy density e in terms of the intensive 
quantities (like T and //) multiplied by their thermodynami- 
cally conjugate extensive densities (like the entropy density s 
and the carrier number density n). Next the equations of mo- 
tion (often conservation laws) for all the extensive densities are 
written down in terms of unknown fluxes and sources. Then 
the volume rate of entropy production S, which is non-zero, 
is written in terms of a divergence and a product of unknown 
fluxes with corresponding gradients of intensive parameters. 
Finally, the condition that 5 > completely determines the 
form of the unknown fluxes (and sources), although the values 
of the transport coefficients are specified only subject to cer- 
tain inequaUties and, for the off-diagonal terms, to the Onsager 
symmetry relations. 

2.1. Rate of Entropy Production 

For a single carrier conductor, we assume that the bulk en- 
ergy density e, bulk number density n, and bulk entropy den- 
sity s are related by 



de = Tds + fidn. 



(4) 



The electrochemical potential jl = fj, — eV, where /i is the 
chemical potential. (In principle, a material can contain multi- 
ple heat-carrying subsystems, which may be at different tem- 
peratures [26], so that de = Tidsi + T2ds2 + . . . ; see Ref. [27] 
for a discussion of the irreversible thermodynamics of such 



©2011 NRC Canada 



Sears, Saslow 



3 



systems.) The continuity relations for the conserved number 
and energy densities are 

fyff 

-Q^+9dr = 0, (5) 
|+5.,f = 0. (6) 

The equation for the non-conserved entropy density is 

ds 

— + dijt =Ss>Q. (7) 

Alternatively one can use the dissipation function 

R = TSs (8) 

as the primary quantity. R has the units of a rate of heating, 
but is only a rate of heating when energy is transformed into 
heat, not when it already is in the form of heat. For a view more 

exclusively based on Sg, see [28]. 

Combining the above equations yields 



(9) 



Following the approach of irreversible thermodynamics, this 
has been written as a single divergence term and the sum of 
products of unknown fluxes with gradients of the thermody- 
namic intensive quantities. (These gradients are also called ther- 
modynamic forces.) 

Since the divergence term may be either positive or negative, 
it must always be zero to ensure that entropy never decreases, 
i.e., Ss > 0. Thus 



and 



< TSs 



(10) 



(11) 



If there is no number current j", then the energy current 
jl consists only of a heat current jf = Tjf. In this case the 
energy is only in the form of heat, and there is no additional 
heat production, although there is entropy production. 

2.2. Bulk Fluxes and Rate of Entropy Production 

In bulk, by the non-negativity of (11), the linearized flux 
densities take the form 



(12) 

(13) 



where k is the thermal conductivity, a is the electrical conduc- 
tivity, e is electric charge, and Lsn = Lns by the Onsager prin- 
ciple. Irreversible thermodynamics cannot provide values for 
any of these material-dependent coefficients, but Kubo theory 
can give these coefficients in terms of equilibrium correlation 
functions [19,20,29]. 



By (1 1)-(13), the rate of entropy production in the bulk (per 
unit volume) is given by 



2^ id.H) (d.T) . (14) 



By 5s > we have k > 0, ct > 0, and L^^ < {aK)/{e'^T). 
As noted in the introduction, pure thermal conduction already 
involves heat flow, so there is no production of heat in that case. 

On the other hand, the entropy production due to current 
flow does cause heating, at the rate per volume of 



L,n{diilMT). 



(15) 



Only the first term is Joule heating. The second term is like 
Thomson heating, in that it can have either sign.^ For a large 
diT, it can even dominate, but the net entropy production Ss 
remains non-negative. 

2.3. Rates of Heating and Entropy Production 

When calculating the rate of heating this elimination, by 
hand, of the part of TSg associated with heat flow is related 
to a similar effect discussed in [30] of damping of a sound 
wave. In that case the mechanical energy E„iech of the sound 
wave (which, implicitly, has zero entropy) dissipates into heat, 
which increases the entropy S of the background system by 
Emech = —TS. Emech is determined by a volume integra- 
tion over the equivalent of R, evaluated for the sound wave, 
and is proportional to the square of the sound wave amplitude 
(including temperature oscillations in the sound wave). This 
results in hot spots (as, for a standing wave, is perhaps famil- 
iar from a microwave oven) that separately diffuse. However, 
this energy is already heat energy. Once deposited as heat, its 
diffusion causes a further increase in entropy, but no additional 
energy goes into the system. Note that in the long-wavelength 
limit the entropy per unit mass ct is conserved, from which one 
can show (see Appendix B) that the entropy of a sound wave 
is zero. 

Another example where entropy increase and heating are 
distinct is a gas of interacting atoms that has a multi-nanometer 
range for repulsion. Let all the atoms initially be placed within 
an interaction volume of one another. When they become ther- 
mally disordered, the increase in entropy can be treated as in 
the present work, but only the interaction energy converts into 
heat. 

In spin-Seebeck experiments [31-33], the applied thermal 
gradient causes heat flow between various subsystems, which 
in turn induces a spin current [27, 34]. Although it increases 
the entropy of the system, the heat flow is not associated with 
heating; however, spin and electrical currents both increase the 
entropy of the system and cause heating. Spin currents are ad- 
dressed in Section 3. 



Because Thomson heating (as discussed in Ref. [11]) involves the 
artificial maintenance of the same temperature distribution both with 
and without current flow, we hesitate to call this cross-term Thomson 
heating, although the latter involves both a temperature gradient and 
a voltage gradient. 
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2.4. Surface Fluxes and Rate of Entropy Production 

At low temperatures, or when the material properties change 
significantly on crossing the interface I, the changes in T and 
ji at the interface are very large, and (14) integrated over the 
surface region (assuming that T and ji are well-defined in this 
region) can be smaller than the surface entropy production rate 
<S/. In this subsection we consider Si. In the next subsection 
we consider the conditions under which Si dominates (14) in- 
tegrated over the surface region. This involves considerations 
of the characteristic mean-free-path / and the distance a over 
which the thermodynamic quantities adjust to the surface. 

The total rate (per unit area) of entropy production at the 
interface. Si, is obtained by integrating the volume rate of en- 
tropy production over the surface region. By (11), taking flow 
only in the x-direction. 



Si = j dxSg = — J 



dx 



(16) 



In steady-state, the energy and number flux densities (10) and 
(13) are uniform across this region. If T and fj, are also nearly 
uniform, by (10) the entropy flux density will also be nearly 
uniform, so 



T 



T 



(17) 



Here (AT) j and (A/i)^^ are the differences of temperature and 
electrochemical potential across the interface region. Johnson 
and Silsbee |6| did not include the effect of variations in the 
chemical potential n, which normally does not matter for met- 
als (which motivated their work). However, variations in ^ of- 
ten are of significance for semiconductors, where they cause 
diffusion. On setting fi, = fi, — eV « —eV, eq. (17) agrees with 
eq. (12) of [6]. See Table 1 for a term-by-term comparison. 

We now apply the same type of irreversible thermodynam- 
ics approach to the surface region as to the bulk. In bulk the 
fluxes are proportional to the gradients of the intensive ther- 
modynamic quantities. At a surface, the fluxes are taken to be 
proportional to the differences across the interface of the inten- 
sive thermodynamic quantities. Thus 



f = -^(AT),-L',„(Am),, 
r = -i;,(AT),-l(AM),. 



(18) 
(19) 



Here, hx is the thermal boundary resistance, and is of the or- 
der of the difference in the products of the specific heat times 
a characteristic sound velocity on each side, and ^ is a sur- 
face conductance, with units l/fl-m^. (We reserve g for the 
g-factor of the charge carriers; for metals g ~ —2.) By the 
Onsager principle (assumed to apply at surfaces as well as in 
bulk), = L'^g. Thus, the total rate of entropy across the 
surface region is 



Si^^ (AT)^ + (A/i)? + 2^ (A/i), (AT), . 



(20) 



The condition Si >0 implies that Hk >0,g>0, and Lf^ < 

{-ghK)/{e'T). 



2.5. Estimates 

We consider a metal-metal interface, for which p, » —eV. 
For characteristic values of current density [35] (J « 10^^ 
A/m^) and surface conductance [5] (g w 10^'' l/$l-m^), a 
characteristic potential difference across the interface is {AV)i 
10-3 V. Then, when (AT)/ = flie appropriate part of TSi 
gives a rate of heating per unit area of 



TZi = TSi = 4(A/i)f « g{AV)j « lO" ^ 



(21) 



We may also apply this to a multilayer system where inter- 
ference effects between the layers, and where mean-free paths 
connecting them, can be neglected [8]. At a separation s of 
100 nm between layers the net effect of the interfaces corre- 
sponds to a bulk conductance of gs w lO^/il-m and the rate 
of heating per unit volume (from the values above) is approxi- 
mately IQi^W/m^. 

As already discussed, there is entropy production because 
of heat flow, but there is no heating rate associated with heat 
flux, because the energy is already in the form of heat. For an 
interface across which there is only a temperature jump, the 
rate of entropy production is given by 



Si 



^(AT)?. 



(22) 



Nevertheless, there is an apparent heating rate, whose value we 
now determine. Typical values for thermal boundary resistance 



[4] (R 



K 



''K 



2 X IQ-^K-m^AV at T = IK for Rh:Fe on 
AI2O3) and energy flux (J^ k. 10-''W/m^) give a value for 
the temperature difference across the interface of (AT)/ w 
2 X lO-^K, so (AT)//T « 0.02. Then 



T5/ = ^(AT)? 



0.2 



W 



(23) 



This apparent heating rate is about ten orders of magnitude 
smaller than for the example of a true surface heating rate due 
to the electrochemical potential gradient, given above in (21). 

2.6. Entropy Production Rates: Bulk vs Surface 

We now consider the conditions under which the interface 
entropy production rate Si can dominate over the near-surface 
space-integral Sb of the bulk entropy production rate Ss- For 
simplicity we consider only carrier flow in the x-direction, with 
no cross-terms, so Lsn ~ 0. 

Three characteristic lengths are associated in this problem: 
one associated with transport (a mean-free path 1), a charac- 
teristic sample size (d), and a distance over which there is an 
interface adjustment (a). 

On equating the bulk and surface carrier currents /uj" and 
using dxjj- ~ (A/i)g /d we find 



(A/2),«-(AA),. 



(24) 



Wifli this result, (15), and (20), for (AT)/ = the integrated 
bulk rate of entropy production near the surface (per unit area) 
is on the order of 



aa- 



( ml 

rf2 



a a 



(25) 
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Thus, when ag/cr 3> 1 (good electrical matching), the contri- 
butions from the bulk electrochemical potential gradients dom- 
inate those near the surface. On the other hand, when ag/a 
1 (poor electrical matching), the contribution from the surface 
jump in electrochemical potential dominates. 

3. Two-Carrier Systems - Spin 

For itinerant magnets the theory should include two carriers. 
A specific case would be the interface between a metal and 
a magnetic material, where spin-up and spin-down electrons 
have different electrochemical potentials on each side of the 
interface. For simplicity we consider that the magnetization di- 
rection M is fixed, and that, for the two adjacent materials, M 
is aligned parallel (Mi, M2 M) or antiparallel (Mi -J> M 
and M2 — M). (Because we do not consider bulk or surface 
spin-transfer torque [36] , we cannot consider non-collinear Mi 
and M2-) Near a magnet, with no external field but with current 
flow, even a nonmagnetic material can develop a nonequihb- 
rium magnetization (also known as spin accumulation). With 
a somewhat different notation, and including terms associated 
with diM, for the bulk many of these results (but not surface 
heating) have been derived previously [36]. Transport of spin 
across surfaces was considered by a number of authors, but 
they did not consider heating rates [6-8]. 

Before introducing the thermodynamics a few definitions 
are needed. First, the theory employs the "magnetization po- 
tential" —H*. H* itself is given by the difference between the 
external fields (magnetic, anisotropy, dipole) and the internal 
field due to exchange. In equilibrium H* = 0? The chemical 
potentials are denoted by ^^4, and are determined, for exam- 
ple, from energy band theory. The electrochemical potentials 
are denoted by fi^^^, and satisfy jx^^^ = /^^^^ — eV. Finally, the 
magnetoelectwchemical potentials are denoted by and 
satisfy 



(26) 



Here, /xq is the permeability of free space, and the gyromag- 
netic ratio 7 = \g\iiB/K where g is the electron spin g-factor 
and /is is the Bohr magneton. In equilibrium jl^ — fii. RecaU 
that [6] does not consider changes in the chemical potential 

'There is some ambiguity in the definition of the external fields and 
the internal fields, but there is no ambiguity in the definition of H* . 
For example, the lattice anisotropy and the dipole fields depends on 
the magnetization, and for that reason can be considered to be in- 
ternal or external. The uniform exchange field is certainly internal, 
but the non-uniform exchange field might be considered internal or 
external. Only the applied magnetic field and the internal exchange 
field should be uniquely considered external and internal. 

''The term magnetization potential is employed by [6] to denote —H*. 
They note that it would be a term in the magneto-electro-chemical 
potential, for which we use 

^Although Johnson and Silsbee [6] formally do not consider changes 
in the chemical potential, in fact such changes are needed very near 
the interfaces, in order that the actual voltage (as opposed to the ex- 
trapolated voltage) be continuous across the interface. This is rele- 



3.1. Rate of Entropy Production 

With these definitions, the bulk energy density, bulk number 
densities for spin up (n^) and spin down (nj,) electrons, and 
bulk entropy density are related by 



de = Tds -\- fi^dn^ -\- fiidni- 



(27) 



Here we neglect a term ^ H ■ dM because we restrict ourselves 
to the longitudinal response. 

The relations for energy density (6) and entropy density (7) 
still apply. Relations for the (non-conserved) number flux den- 
sities are given by 



— +%ti = -^t' 



dt 



(28) 



(29) 



where Sf is the rate at which spin-down electrons flip to spin- 
up electrons. These forms ensure that the total number current. 



is conserved, since 

d (n-f + ni) 



dt 



0. 



However, the dimensionless spin current, 

is not conserved: 

d (n^ - ni) 



dt 



+ diJ[ =2S^. 



(30) 



(31) 



(32) 



(33) 



By equations (6), (7), (28) and (29), we have 



< TSs =di {-3t + Tjt + MtJt. + HJi^) 

- jldiT - - Jlidifil - (/it - P-i) S^. 

(34) 

Again, in order to ensure that the entropy does not decrease, 
the divergence must be zero, so we take 

jt = Tj^ + fifjfi + fiijii. (35) 

Then, 

TSs = -fdiT - j^d,J^ - jidtfii - (fif - pi) Sf. (36) 

Not only does each partial current and the heat current con- 
tribute to the rate of heating, but there is also a contribution 
from spin- flip, involving the difference in magnetoelectrochem- 
ical potential for up and down spins. For conserved spin current 
(i.e., a = 0), eq. (36) agrees with eq. (53) of [6]. See Table 1 
for a term-by-term comparison. 

vant to Fig. 5 of their Appendix. To make the voltage continuous 
at the interface requires charges within a screening length of each 
side of the interface. Each side has its own screening length, and for 
a given side, if the total screening charge (per unit area) is known, 
then the voltage drop near the surface of that side is known. The to- 
tal surface charge per unit area is determined by the difference in the 
electric fields on each side of the interface. The distribution of the 
charges between the two sides is determined by the requirement that 
the voltage be continuous. 
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3.2. Bulk Fluxes and Rate of Entropy Production 

In the bulk, from equation (36), we take the linearized flux 
densities to be 



(37) 
(38) 
(39) 



where and crj, are electrochemical conductivities of up spins 
and down spins, respectively. By the Onsager principle, = 
Lif, Lfs = Lsf, and Lsi — Lis- Further, to ensure the non- 
negativity of (36) even in the absence of gradients of intensive 
variables, Sf is driven by the difference in electrochemical po- 
tentials 



(40) 



where a is proportional to a characteristic spin-flip rate. 
We now expand the magnetoelectrochemical potentials as 



-(0) , 9^^fx - -(0) , '^^■Ia 

Using p,^^ = p,^^^ we then have 
Sf = -a (Ctfef - CiSni) , 



(41) 



(42) 



where Cf = {dfii/dn-i) and Cj. = (dpi/dni) are the inverse 
densities of states. (A similar relation holds in semiconductors, 
where electrons and holes recombine [37].) 

We now relate this to Ti relaxation of the magnetization M 
of a uniform system, where 



M = -M^(n,-n,). 



(43) 



Recall that g is the charge carrier g-factor. Particle conservation 
gives Sni = —Sn^f, so 5M = — l^l/iB^n-]- and = —a{C-f + 
Ci)Sn-^. Then (33), with 9, neglected, yields 



dM , , d , M 
^ = -|5K^(<5n,) = --, 



(44) 



where Tf^ = a{C-f + Ci). 

The dissipation function R = TSs is given for the bulk by 



TSs =^ {d,Tf + ^ (d^Ptf + ^ {d,Pi) 

+ 2L,t i^^nt) {d,T) + 2Lsi (a,/2j (S^T) (45) 

+ 2iti {dtfH) {dtp,i) + a (/i^ - f^ f . 

This is sufficiently complex that each term deserves comment. 

The term in (diT)"^ is from heat current, the terms in {dipf)"^ 
and (dipi)'^ are from Joule losses of the individual carriers 
[38], the next three are cross-terms, and the last term gives the 
dissipation due to spin-flip processes. Ss > forces various 
conditions on both the diagonal and the cross-terms, of which 
the latter usuaUy are smaU. The diagonal terms satisfy k > 0, 
> 0, (Tj. > 0, and a > 0. Note that it is not the current 



or spin current (both of them thermodynamic fluxes) that de- 
termines the rate of entropy production and heating, but rather 
the gradients of the magnetoelectrochemical potentials (both 
of them thermodynamic forces). 

Ref. [39] considers bulk entropy production due to spin ac- 
cumulation in a non-magnetic material that is adjacent to a fer- 
romagnet. The present work considers the bulk entropy pro- 
duction in both the non-magnetic material and the ferromag- 
net, as well as the total rate of entropy production at and near 
the interface (see below). 

3.3. Dissipationless Spin Currents 

In recent years there has been considerable interest in the 
possibility of spin currents that do not cause dissipation. Al- 
though in [40,41] an electric field is applied, Rashba notes that 
it is possible to have (equilibrium) spin currents, even without 
an applied electric field [42] . This situation can be incorporated 
within the present theory by: 



(a) letting j^i jjf ^ + Sj-^, and jj,, j|f ^ + Sj^, in 
(28) and (29) and the remainder of that subsection, with 
0; and 



d4f 



(b) taking (38) to apply to and (39) to apply to ^in- 

Because the rate of entropy production involves the thermo- 
dynamic forces, (45) still applies. For the rate of heating, one 
should apply (36) with the jfdiT term omitted. 

3.4. Surface Rate of Entropy Production 

Consider flow along x, so we may drop the directional in- 
dices on the fluxes. Then, rewriting (36) with (30) and (32) 
gives 



(46) 



In steady-state J" is constant across the interface region, but 
is not. Moreover, is not obviously a near-conserved quan- 
tity, unlike for a single carrier. However, since is conserved, 
we use (35) to write 



T 2 \T t) 2 VT T 



(47) 



With both J" and the difference in magnetoelectrochemical 
potentials considered to be first order in smaU deviations from 
equilibrium, the last term is second order. Thus actuaUy is a 
near-conserved quantity. 

Integrating the volume rate of heating (46) over the interface 
region yields 

TSi = - f (AT), i J" ((Amt), + (A/x;),) 

(48) 
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Integrating the third term by parts and using equation (33) with 
the time-derivative set to zero (steady-state, so that dxJ" = 
25) gives 

~\ [ dxJ''d:,{ji^-Jl^) 

(49) 

= j^dx (Mt - H) - \^ {J" (At - H))i ■ 

The first term on the right-hand-side (RHS) of (49) cancels the 
last term on the RHS of (48), and the second term is evaluated 
on each side of the interface. Then (48) becomes 

TSi = - f (AT), - I J" ((A/i^), + (A/2J,) 

1 ^ (50) 

- 2^ ("^" (At - ■ 

The last term appears to be well-defined, but because J"^ is not 
conserved, it is not clear how to interpret this unambiguously. 
Nevertheless, if the spin diffusion length, over which up and 
down spins flip, is sufficiently long relative to the surface re- 
gion, so that one can measure J"^ within a spin diffusion length 
of the interface, then this last term should be clearly defined 
(otherwise it may not be measurable).^ In what follows we will 
assume that this holds. For conserved spin current, eq. (50) 
agrees with eq. (44) of [6]. See Table 1 for a term-by-term 
comparison. 

3.5. Surface Fluxes and Rate of Entropy Production 

In the spirit of irreversible thermodynamics, the entropy and 
the spin up and spin down number fluxes can be hnearized in 
differences in the appropriate intensive thermodynamic quan- 
tities across the interface, so 

f = -^(AT),-L',^(A/2t),-i',;(A/24),, (51) 

it = -4,, (AT), - ^ (AMt), - (AMi), , (52) 
h = -^1. (AT), - Lit (A/i^), - |i (Am^), . (53) 

Here and are surface conductances of spin up and spin 
down particles, and by the Onsager principle L!j.j^ = L'^^,L'^g = 
L^t' = ^'is- ^ calculation of a spin-dependent in- 

terfacial surface resistance, see Ref. [43]. 

From (30) and (32) the total number current and spin current 
can be written as 

J" = - + L[,) (AT), - + L[,) (A/lt), 

-(|+^U)(AAJi' (54) 

^Because /I^.J. vary exponentially within each material's spin- 

diffusion length from the interface, there are four intervals over 
which TSi of (50) is clearly defined. Taking the edges of the inter- 
vals to be either much closer or much farther than the spin-diffusion 
lengths from the interface yields four possible intervals. 



J- = - (L^, - L[,) {AT), - (I - L[,) (A/2t), 

+ (|-^n)(AA^)/- (55) 
Substitution of the currents in (51) and (54) into (50) yields 

TSi = -^A (r (Mt - Mi)), + ^ (AT)? 

+ ^(|H-^lt)(AMt)? + ^(|+ilt)(AM.)? 
+ ^(34t + ig (AMt),(AT), 

+ i(3L:,, + Tg (A^,), (AT), 

+ 1 + 2i^lt) (AAt)/ (Ami), • (56) 

If the spin current J'^ is approximately uniform near the sur- 
face, then use of (55) gives that (56) simplifies to 

TSj =^ (AT)? + g (A^t)? + (A/2;)? 

+ 2L;t (A^t)/ (AT), + 2L;^ (A/xJ, (AT), 

+ 2Lit(AMt)/(A/ii),. (57) 

as in [6], which gives an approximation for each of the coeffi- 
cients. 

3.6. Comparison of Electric and Spin Current Heating 

Eq. (57) permits a comparison of surface heating due to elec- 
tric current (equivalently, due to a voltage jump across the sur- 
face region) with heating due to spin current (equivalently, due 
to the difference in //oiT* • M across the surface region). We 
consider a metal-metal interface where temperature is uniform, 
spin is conserved across the interface, and fif-f ^ 54, ~ 5/2. 

For the purposes of estimation, we neglect the chemical po- 
tentials ji^ and ji]^. The limitations of this approximation are 
discussed above. Then, (26) gives 

(AMt,i)?--e^(Ay)? + :i^Mg(AiJ*)? 

Tineiia{AV)i{AH*)i. (58) 

Here we define H* = H* ■ M . To find the heating due to 
(Ay)/ and /xo(A-ff*)/, we substitute (58) into the second and 
third terms of (57). 

Suppose (as above) that the surface region is characterized 
by surface conductivity g = 10^'' l/fl-nP, and voltage dif- 
ference {AV)i = 10~^ V. Neglecting cross-terms, the rate 
of heating, per unit area, due only to the voltage difference is 
given by 

j^elec ^ j,gelec ^ ^ [{AV)if « 10^ (59) 

(Cancellation of the term proportional to (AF)/(AiI*)/ is 
only approximate.) 
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On the other hand, neglecting cross-terms, the rate of heat- 
ing due only to spin current is given by 



^72 



Jospin rj-i 

=^[\g\t^BHo{AH*)i] 
For 1^1 w 2 and /is » 5.8 * 10"^ eV/T, 
Tlf ™ « 3.4 * 10*^ [no{AH*y^ ^ 



T2-] 



n2 • 



(60) 



(61) 



Thus, a iJ,o{AH*)i « 20 T gives about the same heating as a 
voltage difference of (AV) / = IQ-^ V. 

4. Summary and Conclusions 

We have applied the methods of irreversible thermodynam- 
ics to study the rate of entropy and heat production in bulk 
and associated with an interface, first for an ordinary conduc- 
tor (this includes insulators) and then for a conducting magnet 
whose magnetization direction M is fixed. In addition to mag- 
netic metals, the present results apply to magnetic insulators 
and semiconductors due to the incorporation of the chemical 
potential, neglected previously. We also show how equilibrium 
spin currents can be included in the present theory. This is rel- 
evant to heating that takes place in nanoelectronic systems and 
spintronic systems, especially in multilayers. 

We conclude by noting that the present work is related to 
recent work on "spin caloritronics," whereby heat currents can 
cause spin currents and spin currents can cause heat currents 
[44^6]. The former, known as the spin-Seebeck effect, has 
recently been measured [31-33] using the inverse spin Hall 
effect, and in one case [33] displays a profile that is associated 
with spatially exponential decay away from the heat input and 
output leads [27,34]. 
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Appendix A: Entropy-Maximum-Based 
Variables 

One may write the increase in entropy in terms of entropy- 
maximum-based intensive variables (e.g., 1/T, -/i/T). This is 
most easily shown by writing the time derivative of the entropy 
density as 



Eds . 



(62) 



where the index k here is used to denote densities Xk of the 
entropy-maximum-based extensive variables. These include the 
energy density e, which is conjugate to 1/T. We take an equa- 
tion of motion that includes both flux jf and source Sk (Sk = 
if Xfe is conserved, but for spin-flip scattering it is important): 



Xk = -diJi +Sk, 
Then (62) gives 



dxk 



(63) 



(64) 



Solving for Ss (which must be non-negative) and rewriting so 
that there is an independent gradient term yields 



0<Ss =di 



J! 



ds 



' dxk 



(65) 



Since Ss > Q but the divergence may take either sign, the 
argument of the divergence must be zero (modulo a curl, which 
we neglect). This relates the entropy flux to the other fluxes, but 
does not determine them. 
To obtain the other fluxes, we set the divergence to zero: 



ds 
dxk' 



(66) 



For sound in the long-wavelength limit a is constant, so da = 
0. Moreover, Sp = 0. Therefore the first four terms are zero. 
Further, the entropy is determined from the rest frame, and thus 
has no velocity-dependence, so the sixth term is also zero. Fi- 
nally, since s = pa and each derivative is taken with the other 
variables held constant, the fifth term also is zero: 



dp^ 



da 
'dp 



= 0. 



(68) 



Therefore, to second order in fluctuations, a sound wave has 
zero entropy; its energy is solely in the form of potential energy 
and kinetic energy. 

Appendix C: Notationai Variation 

CaUen (see eq. (17.3) of Ref. [11]) uses gradients of en- 
tropic variables (e.g., 1/T and —ji/T) for generalized thermo- 
dynamic forces. Except for Appendix A, the present work uses 
gradients of energy-minimum variables (e.g., T, p). Morse (see 
eqs. (8.36) and (8.48) of Ref. [12]) shows the results of both 
approaches. Johnson and Silsbee (see eq. (12) of Ref. [6]) use 
neither convention, and have gradients of 1/T and voltage V 
for generalized forces. 

These differences make comparison complex. Table 1 com- 
pares the quantities of the present work to those of Ref. [6]. 
These authors take the electron charge ase = — |e|,^ = hb 
as the Bohr magneton, and the spin g-factor as positive. We 
take the gyromagnetic ratio to be7 = \g\pB/^- Recall that 
Ref. |6| neglects the chemical potentials {p, pf and pi) rel- 
ative to —eV. Thus Table 1 applies for p. « —eV. 



This is the entropy-maximum variable analog of (1 1) and (36), 
which are in energy-minimum variables. Only the first term 
would appear if the extensive quantities were conserved; the 
second term represents the effect of sources. To obtain the irre- 
versible thermodynamics one can expand the unknown sources 
and fluxes in terms of the presumably known ds/dxk's and 
their gradients. 



Appendix B: Sound Waves 

Consider a sound wave. Let subscript denote an equilib- 
rium value. With a = s/p, the time-averaged change in en- 
tropy density due to a sound wave is 

ds-^ ds-^ 1 d'^s . , 1 d^s ,„ 

1 ds 
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Table 1. The relationship between quantities in the present work and in Ref. [6]. We take e > in this table and in the present work, 
whereas Ref. [6] takes e < 0. 





Bulk (Continuous System) Quantities 


Surface (Discrete System) Quantities 




Present 


Ref. [6] 


Present 


Ref. [6] 


Entropy Production 




Suniv 


Si 




Rntronv Pliix 


f 


Jq 
T 


f 


T 


Charge-Carrier Flux 


r, J" 


J 

-11. 
e 


r, J" 


r 

-h. 

e 


Spin Flux 


J" 


2Jm 




27m 

7ft 


i hermal Conductivity 
and Conductance 


K 


^22 
y2 




^22 

rp2 


Electrical Conductivity 
and Conductance 


a 


Ln 


S 


-'^11 

T 


Electrical Conductivity and 


(7f 


ill e^Lss 61/31 
4T 72/|2r --ihT 


5t 


i'n e^L'ys CL31 
4r 72^27 7ftT 


Electrical Conductivity and 
fntiHiiPtfinpf* fnr 1 -SlniTis 




ill e^Laa eLai 




i'll e2L33 61/31 

4T 72/i2T 7ST 


son Onsager Coefficient 


Lsti 


L21 


r' 


r' 

-1^21 

eT2 


s -o^ t Onsager Coefficient 




L32 1 il2 

7fiT2 2eT2 


Tl 


7"' T' 
^32 1 ^12 

7fiT2 2eT2 


s O 4. Onsager Coefficient 


Lai. 


L32 1 il2 

7/IT2 2eT2 


i'si 


r' T' 

^32 . ^12 

7/iT2 2eT2 


t 4- Onsager Coefficient 




ill L33 




in iaa 


4e2T 72ft2T 


4e2T 72^27 
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